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ABSTRACT 

Bulk dark energy properties are determined by the redshift evolution of its pressure-to-density ratio, 
w^c{z). An experimental goal is to decide if the dark energy is dynamical, as in the quintessence (and 
phantom) models treated here. We show that a three-parameter approximation w^e^z] Eg, £4,00, Cs) fits 
well the ensemble of trajectories for a wide class of late-inflaton potentials V{4>). Markov Chain Monte 
Carlo probability calculations are used to confront our Wdeiz) trajectories with current observational 
information on Type la supernova. Cosmic Microwave Background, galaxy power spectra, weak lensing 
and the Lyman-a forest. We find the best constrained parameter is a low redshift slope parameter, 
Es oc {d In V/dcf))"^ when the dark energy and matter have equal energy densities. A tracking parameter 
£000 defining the high-redshift attractor of 1 -l-Wde is marginally constrained. Poorly determined is Cs, 
characterizing the evolution of Eg , and a measure of d^ In V/dcjP' . The constraints we find already rule 
out some popular quintessence and phantom models, or restrict their potential parameters. We also 
forecast how the next generation of cosmological observations improve the constraints: by a factor of 
about five on e^ and e^oo, but with Cs remaining unconstrained (unless the true model significantly 
deviates from ACDM). Thus potential reconstruction beyond an overall height and a gradient is not 
feasible for the large space of late-inflaton models considered here. 
Subject headings: cosmology, dark energy, quintessence 



1. INTRODUCTION 
1.1. Running Dark Energy and its Equation of State 

One of the greatest mysteries in physics is the nature 
of dark energy (DE) which drives the present-day cos- 
mic acceleration, inferred to exist from su pernovae data 
(|Riess et all 119981 : iPerlm utter et al."1999'). and from a 
combination of cosmic microwave backgrou nd and large 
scale structure data (JBond fc Jaffd 119991 ). Although 
there have been voluminous outpourings on possible the- 
or etical explanations of d ark energy, recently r eviewed 
in'Copeland et al.' C2006') : 'Kamionkowskil ()2007[) : iLindeii 
dlOOSa^b); Silvcstri & Trodden (200^rwe are far from 
consensus. An observational target is to determine if 
there are temporal (and spatial) variations beyond the 
simple constant A. Limits from the evolving data con- 
tinue to roughly center on the cosmological constant case 
A, which could be significantly strengthened in near- 
future experiments - or ruled out. We explore the class 
of effective scalar field models for dark energy evolution 
in this paper, develop a 3-parameter expression which 
accurately approximates the dynamical histories in most 
of those models, and determine current constraints and 
forecast future ones on these parameters. 

The mean dark energy density is a fraction ilde of the 
mean total energy density, 

Pdo(lna) EE ptot^dc , 3A/p^i7^(a) = ptot(a) , (1) 

which is itself related to the Hubble parameter H through 
the energy constraint equation of gravity theory as in- 
dicated. fide(o) rises from a small fraction relative to 
the matter at high redshift to its current ~ 0.7 value. 
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Here Afp = l/V8iG^ = 2.44 x lO^^GeV, is the reduced 
Planck mass, Gm is Newton's gravitational constant and 
h and c are set to unity, as they are throughout this 
paper. 

Much observational effort is being unleashed to deter- 
mine as much as we can about the change of the tra- 
jectory pdc with expansion factor a, expressed through a 
logarithmic running with respect to In a: 



Idlnpdc ., ^ 3 

o-7j = edc Ina) = - 1 
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This pdc-i'un is interpreted as defining a phenomenolog- 
ical average pressure-to-density ratio, the dark energy 
equation of state (EOS). The total "acceleration factor", 
e = 1 + q, where the conventional "deceleration param- 
eter" is g = —aaja} = —dhi{Ha)/dliia, is similarly re- 
lated to the running of the total energy density: 
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- edo^dc , ^m + f^dc — 1 • (3) 



Eq. ^ shows e is the density-weighted sum of the ac- 
celeration factors for matter and DE. Matter here is ev- 
erything but the dark energy. Its EOS has e™ = 3/2 in 
the non-relativistic-matter-dominated phase (dark mat- 
ter and baryons) and e^ = 2 in the radiation-dominated 
phase. 

1.2. The Semi-Blind Trajectory Approach 

We would like to use the data to constrain pdc(lna) 
with as few prior assumptions on the nature of the tra- 
jectories as is feasible. However, such blind analyses are 
actually never truly blind, since pdc or Wdc is expanded 
in a truncated basis of mode functions of one sort or 
another: necessarily there will be assumptions made on 
the smoothness associated with the order of the mode 
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Fig. 1. — The marginalized 68.3% CL (inner contour) and 95.4% 
CL (outer contour) constraints on wq and «)a for the conventional 
DETF parainetrization w = wq + Wa{l — a), using the current 
data sets described in § |4] The white point is the cosmological 
constant mod el. The soli d red line is a slow-roll consistency rela- 
tion, Eq. I|50| l derived in 5 16.21 (for a fixed QmO = 0.29, as inferred 
by all of the current data). The tilted dashed gray line shows 
Wa = — 1 — Wo- Pure quintessence models restrict the parameter 
space to 1 -|- lUQ > and Wa above the line, whereas in the pure 
phantom regime, 1 + too < ~0 and Wa would have to lie below the 
line. Allowing for equal a priori probabilities to populate the other 
regions is quite problematical theoretically, and indeed the most 
reasonable theory prior would allow only the pure quintessence do- 



expansion and on the measure, i.e., prior probability, of 
the unknown coefficients multiplying the modes. The 
most relevant current data for constraining Wdo, SNIa 
compilations, extends back only about one e-folding in 
a, and probes a double integral of Wdc, smoothing over 
irregularities. The consequence is that unless Wde was 
very wildly varying, only a few parameters are likely to 
be extractable no matter what expansion is made. 

Low order expansions include the oft-used but non- 
physical cases of cons tant wq j^ —1 and the 2-parameter 
linea r expansion in a (ICh cvallicr & Polarski 2001; Lindcr 
200l IWrng fc Tegma^k 2004- Upadhve et al. 2005: 



ZW^ IWang &: legmarK: zom- Upadlrve et al. ZOOb; 
Liddle et al.ll2006HFrancis et al.ll2007tlWang et al.ll2007l) . 



w{a) = wq +Wa{l - a) , 



(4) 



adopted by the Dar k Energy Task Force (DETF) (see 
lAlbrecht et al.l 120061) . The current observational data 
we use in this paper to constrain our more physically- 
motivated parameterized trajectories are applied in Fig- 
urc[l]to Wq and Wa, assuming uniform uncorrelated priors 
on each. The area of the nearly-elliptical 1-sigma error 
contour has been used to compare how current and pro- 
posed dark energy probes do relative to each other; its 
inverse defines the DETF "figure of merit" (FOM). 

Constant Wde and Eq. ^ can be considered to be 
the zeroth and first order polynomial expa nsion in the 
variable 1 — a. Why not redshift z (Up adhve et al.l 
[2005; Lindcr & Hutcrcr 200j|) or the scale factor to some 



Jloi 



power (|Linder fc Huterer.,20051 ) or In a? Why pivot about 
a = l. 

Why only linear in a? Why not expand to higher or- 
der? Why not use localized spline mode functions, the 



simplest of which is a set of contiguous top hats (unity 
in a redshift band, zero outside)? These cases too have 
been much explored in the literature, both for current 
and forecast data. The method of principal component 
(parameter eigenmode) analysis takes orthogonal linear 
combinations of the mode functions and of their coeffi- 
cient s, rank-ordered by the error matrix eigenmodes (see 
e.g.. iHuterer fc StarkmanI 120031 : iCrittenden fc PogosianI 
120051) ^ As expected only a few linear combinations of 
mode function coefficients are reasonably determined. 
Thus another few-parameter approach expands Wdc in 
many modes localized in redshift but uses only the low- 
est error linear combinations. At least at the linearized 
level, the few new parameters so introduced are uncorre- 
lated. However, the eigenmodes are sensitive to the data 
chosen and the prior probabilities of the mode parame- 
ters. An alternative to the w^c expansion is a direct mode 
expan sion in e, considered at low order bv lRapetti et all 
()2007|) : since the Cm part is known, this is an expansion 
in (edo — em)^de, in which case Wde becomes a derived 
trajectory; for Wde parameterizations it is e(a) that is the 
derived trajectory. 

The partially blind expansions of 1 -f Wdc are similar 
to the early universe infiation expansions o f the s calar 
power spectrum Vsi^nk), e.g., iBond et all (|2010[ ). ex- 
cept in that case one is "running" in "resolution". Ink. 
However there is an approximate relation between the 
comoving wavenumber k and the time when that specific 
fc-wave went out of causal contact by crossing the instan- 
taneous comoving "horizon" parameterized by Ha. This 
allows one to translate the power spectra In fc-trajectories 
into dynamical Inifa-trajectories, Vs(}nHa). Lowest or- 
der is a uniform slope n^ — 1 = dlnPs/ In fc; next order is 
a running of that slope, rather like the DETF linear Wde 
expansion. Beyond that, one can go to higher order, e.g., 
by expanding in Chebyshev polynomials in In fc, or by us- 
ing localized modes to determine the power spectrum in 
fc-bands. Generally there may be tensor and isocurva- 
ture power spectra contributing to the signals, and these 
would have their own expansions. 

In inflation theory, the tensor and scalar spectra are 
related to each other by an approximate consistency con- 
dition: both are derivable from parameterized accelera- 
tion factor trajectories defining the inflaton equation of 
state, 1 -|- 'w^{a) = 2e/3 (with zero $7^ in early universe 
inflation). Thus there is a very close analogy between 
phcnomenological treatments of early and late universe 
inflation. However, there is a big difference: power spec- 
tra can be determined by over ~10 e-foldings in k from 
CMB and Large Scale Structure data on clustering, hence 
over ~10 e-foldings in a, whereas dark energy data probe 
little more than an e-folding of a. This means that higher 
order partially blind mode expansions are less likely to 
bear fruit in the late-inflaton case. 

1.3. Physically-motivated Late-inflaton Trajectories 

This arbitrariness in semi-blind expansions motivates 
our quest to find a theoretically-motivated prior prob- 
ability on general dark energy trajectories character- 
ized by a few parameters, with the smoothness that is 
imposed following from physical models of their origin 
rather than from arbitrary restrictions on blind paths. 
Such a prior has not been much emphasized in dark en- 
ergy physics, except on an individual-model basis. The 



theory space of models whose trajectories we might wish 
to range over include: 

(1) quintessence, with the acceleration driven by a scalar 
field minimally coupled to gravity with an effective 
potential with small effective mass and a canonical 
kinetic energy (iRatra fc Peebles! 119881 : iWette rich 1988; 



Frieman et al 
2000[ 



19951: 



Binctruv 1999; Barreiro et al 
1999; _ .Copcland et ali ,2000: 



IBrax fc Martini , 

de La Macorra fc Stephan-Ottol 120011 : iLindei 



20061: 



Huterer fc Peirisll2007l : ILinderil2n08l ) 

(2) k-essence, late inflation with a non-canonical 
scalar field kinetic energy a s well as a po t ential 
(lArmendariz-Picon et alJ 120011: iMalauarti et~all 120031: 

iMalauarti et al.ll2003t Ide Putter fc Lindeiil2007l) : 

(3) F{R, (j)) models, with Lagrangians involving a 
function of the curvature scalar _R or a general- 
ized dilaton 0, that can be transformed through 
a conformal mapping into a theory with an effec- 
tive scalar field, albeit with interesting fifth force 
couplings to the matte r sector poppi ng out in the 
transform (ICapozziello"fc Fang 2002; Car roll et aLll2004l: 



de La Cruz-Dombriz fc Dobado 2006; Noiiri & OdintsoyI 
2006; dc Felice fc Tsuiikawa 2010; Sotiriou & Faraoni 

(4) phantom energy with Wde < ~1, an d negative effec 



five kinetic energy, also with a potential (iCopeland et"~al 



mm 



20061 ICarrohet al1[200l ICaldwelll [20021: ICaldwell eTa! 



In this paper, we concentrate on the quintessence class 
of models, but allow trajectories that would arise from 
a huge range of potentials, with our fitting formula only 
deviating significantly when |1 -I- Wdo| > 0.5. A class of 
models that we do not try to include are ones with dark 
energy having damped oscillations before settling into a 
minimum A (Kaloper & Sorbo 2006; Dutta & Schcrrer 
[20081) . 

We also extend our paths to include ones with Wdc < 
— 1 , the phantom re gime where the null energy condition 
([Carroll et al.l I2OO30 is violated. This can be done by 
making the kinetic energy negative, at least over some 
range, but such models are ill-defined. There have been 
models proposed utilizing e xtra UV physics, such as with 
a Lorentz- violating cut-off (iCline et al.l[2004D or with ex- 
tra fields ([Sergienko fc Rubakovl |2008[ )^ Regions with 
Wdc < ^1 are consistent with current observations, so 
a semi-blind phenomenology should embrace that pos- 
sibility. It is straightforward to extend our quintessence 
fitting formula to the phantom region, so we do, but show 
results with and without imposing the theoretical prior 
that these are highly unlikely. 

The evolution of Wdc for a quintessence field (f) de- 
pends on its potential V{4i) and its initial conditions. 
For a flat Friedmann-Robertson- Walker (FRW) universe 
with given present matter density pmo and dark en- 
ergy density pdoO: any function u'de(a) satisfying < 
1 -|- Wde ^ 2 can be reverse-engineered into a poten- 
tial V(0) and an initial field momentum, Il^^ini = ipini- 
(The initial field value (/)ini can be eliminated by a 
translation in V{(l)), and the sign of Il^^ini by a re- 
flection in l/((/i).) Given this one-to-one mapping be- 
tween {V{(j)),\ij>ini\;p„ia} and {w(a), Pdoo; Pmo}, we may 
think we should allow for all Wdc trajectories. Gener- 
ally, these will lead to very baroque potentials indeed. 



and unrealistic initial momenta. The quintessence mod- 
els with specific potentials that have been proposed in 
the literature are characterized by a few parameters, 
and have a long period in which trajectories relax to at- 
tractor ones, leading to smooth trajectories in the ob- 
servable range at low redshift. These are the simple 
class of potentials we are interested in here: our phi- 
losophy is to consider "simple potentials" rather than 
simple mathematical forms of Wde(a)- The potentials 
which we qua ntitatively discuss in later sections include: 
powe r-laws ([Ratra fc Peebles 1988; S tcinh ardt et al] 
Il999| ): expon entials ([Ratra fc Peebles! 1 1 9881): double 
exponen tials (iBarreiro et al! [20001: [ Sen fc Sethil [20021: 
[Neupan e 2004allb!: [Jarv et al.[[2004!) : cosines from pseudo- 
Nambu Goldstone bosons (pnGB) ([Frieman et al.![r995[ : 
Kaloper fc Sorbo 2006. ) ; supergravity (SUGRA) mo- 
tivated models (Brax fc Martini Il999f ) . Our general 
method applies to a much broader class than these 
though: we find the specific global form of the potential 
is not relevant because the relatively small motions on its 
surface in the observable range allow only minimal local 
shape characteristics to be determined, and these we can 
encode with the three key physical variables parameter- 
izing our Wdc- Even with the forecasts data to come, we 
can refine our determinations but not gain significantly 
more information - unless the entire framework is shown 
to be at variance with the data. 

1.4. Tracking and Thawing Models 

Quintessence models are often classifie d into track- 
ing and thawing models fCaldwell fc Lind eill2005l: iLindeii 
20061 [2008t IHuterer fc Peiria,2007; .Sullivan et al.! 120071: 
Scherrer fc Senll2008D . Tracking models were first pro- 
posed to solve the coincidence problem, i.e., why the dark 
energy starts to dominate not very long after the epoch 
of matter-radiation equality ('^ 7.5 e- folds). However, a 
simple negative power-law tracking model that solves the 
coincidence problem predicts Wdc j^ ~0.5 today, which is 
strongly disfavoured by current observational data. In 
order to achieve the observationally favoured Wdc ~ — 1 
today, one has to assume the potential changes at an en- 
ergy scale close to the energy density at matter-radiation 
equality. The coincidence problem is hence converted to 
a fine-tuning problem. In this paper we take the track- 
ing models as the high-redshift limit, and parameterize 
it with a "tracking parameter" e^oo that characterizes 
the attractor solution. The thawing models, where the 
scalar field is frozen at high redshift due to large Hub- 
ble friction, can be regarded as a special case of tracking 
models where the tracking parameter is zero. Because 
of the nature of tracking behaviour, where all solutions 
regardless of initial conditions approach an attractor, we 
do not need an extra parameter to parameterize the ini- 
tial field momentum. 

In either tracking models or thawing models the scalar 
field has to be slowly-rolling or moderately-rolling at low 
redshift, so that the late-universe acceleration can be 
achieved. The main physical quantity that affects the 
dynamics of quintessence in the late-universe accelerat- 
ing epoch is the slope of the potential. The field momen- 
tum is "damped" by Hubble friction, and the slope of the 
potential will determine how fast the field will be rolling 
down; actually, as we show here, it is the ^-gradient of 
In V that matters, and this defines our second "slope pa- 



rameter" Ss- 

It turns out that a two-parameter formula utilizing just 
£0oo and Eg works very well, because the field rolls slowly 
at low redshift, indeed is almost frozen, so the trajec- 
tory does not explore changes in the slope of In V. How- 
ever, we want to extend our formula to cover the space 
of moderate-roll as well as slow-roll paths. Moreover, 
in cases in which In V changes significantly, even slow- 
roll paths may explore the curvature of In V as well as 
its (/)-gradient. Accordingly, we expanded our formula 
to encompass such cases, introducing a third parameter 
Cs, which is related to the second 0-derivative of InT^, 
explicitly so in thawing models. 

1.5. Parameter Priors for Tracking and Thawing 
Models 

Of course, there is one more important fourth param- 
eter, characterizing the late-inflaton energy scale, such 
as the current dark energy density pdc,o at redshift 0. 
This is related to the current Hubble parameter Hq = 
100ft. km s~^Mpc~"'^, the present-day fractional matter 
density ^mo j and the current fractional dark energy den- 
sity r^A = ridcfi = 1 - flmoh'^/h^: Pdc,o = SMp'^H^riA- 
With CMB data, the sound-crossing angular scale at re- 
combination is actually used, as is the physical density 
parameter of the matter, ftmoh'^, which is typically de- 
rived from separately determined dark matter and bary- 
onic physical density parameters, fldmo^^ and ftboh^- 
With any nontrivial lUdc EOS, this depends upon the pa- 
rameters of Wdc as well as on one of h or the inter-related 
zero redshift density parameters. The important point 
here is that the prior measure is uniform in the sound- 
crossing angular scale at recombination and in ildmo^^ 
and ftboh^, not in h or the energy scale of late-inflation. 
If the quantities are well-measured, as they are, the spe- 
cific prior does not matter very much. 

What to use for the prior probability of the parameters 
of Wdo, in which most are not well-determined? The prior 
chosen does matter, and has to be understood when as- 
sessing the meaning of the derived constraints. Usually 
a uniform prior in wq and Wa is used for the DETF form, 
Eq. ^. The analogue for our parameterization is to be 
uniform in the relatively well-determined £=,, but we show 
what variations in the measure on it do, e.g., using y/\e^ 
or disallowing phantom trajectories, Ss > 0. The mea- 
sure for each of the other two parameters is also taken to 
be uniform. (Our physics-based parameterization leads 
to an approximate "consistency relation" between wq and 
Wq, a strong prior relative to the usual uniform one.) 

In § [21 we manipulate the dynamic equations for 
quintessence and phantom fields to derive our approx- 
imate solution, Wdc{a\£s,£ci>oc,Cs,^mo)- We show how 
well this formula works in reproducing full trajectories 
computed for a variety of potentials. In § 2] we de- 
scribe our extension of the CosmoMC program package 
(jLewis fc B ridle 200^ to treat our parameterized DE 
and the following updated cosmological data sets: Type 
la Supernova (SN), galaxy power spectra that probe large 
scale structure (LSS), weak lensing (WL), Cosmic Mi- 
crowave Background (CMB), and Lyman-a forest (Lya). 
The present-day data are used to constrain our dynam- 
ical Wdcici). In § [5l we investigate how future observa- 
tional cosmology surveys can sharpen the constraints on 



Wdcia\es,Scj,oo,Cs,^mo)- We restrict our attention to fiat 
universes. We discuss our results in § |6l 

2. LATE-INFLATION TRAJECTORIES AND THEIR 
PARAMETERIZATION 

2.1. The Field Equations in Terms of Equations of 

State 

We assume the dark energy is the energy density of a 
quintessence (or phantom) field, with a canonical kinetic 
energy and an effective potential energy V{4>). This is 
derived from a Lagrangian density 



£ = ±ld,^d>^<t> - vi4>) , 



(5) 



where -I- is for quintessence (— for phantom). Through- 
out this paper we adopt the {+,—,—,—) metric signa- 
ture. For homogeneous fields the energy density, pressure 
and equation of state are 



P^ = \nl + v{cj,), p^^^uyv {<!>), 


(6) 


_P^+P^ n2 d<f> 


(7) 


PS Pd, at 



Since we identify the dark energy with an infiaton, we 
hereafter use (j> as the subscript rather than "de". The 
field (j) does not have to be a fundamental scalar, it could 
be an effective field, an order parameter associated with 
some collective combination of fields. A simple way to 
include the phantom field case is to change the kinetic 
energy sign to minus, although thereby making it a ghost 
field with unpalatable properties even if it is the most 
straightforward way to get 1 -|- Wde < 0. 

The two scalar field equations, for the evolution of 
the field, (j) — H^, and of the field momentum, H^ = 
SHU^ — dV/d(j), transform to 






id-^/dlna , i}} = (/)/(V2Mp) 
dip 




(8) 



(9) 



The sign is arbitrary, depending upon whether the field 
rolls down the potential to larger ip {+) or smaller ip (— ) 
as the universe expands. For definiteness we take it to be 
positive. The second equation can be recast into a first 
order equation for d In U e^^^j d In a which is implicitly 
used in what follows, we prefer to work instead with a 
first integral, the energy conservation equation ([2]). 

As long as H^ is strictly positive, tp is as viable a time 
variable as In a, although it only changes by y^^H^ in 
an e-folding of a. Thus, along with trajectories in In a, 
we could reconstruct the late-inflaton potential V{<p) as 
a function of (p and the energy density as a function of 
the held: 



(11) 
(12) 




These "constraint" equations require knowledge of both 
60 and 57^, but the latter runs according to 



2 din a 



(13) 



hence is functionally determined. Eq. (J13p is obtained by 
taking the difference of Eqs. ^ and ^ and grouping all 
Q^ terms on the left hand side. 

For early universe inflation with a single inflaton, D,^ = 
1, and £0 = ££ = e. The field momentum therefore 
obeys the relation II^ = —2Mp^dH/d(j}, which is derived 
more generally from the momentum constra int equation 
of general relativity (jSalopek fc Bondlll99C0 . 

2.2. A Re-expressed Equation Hierarchy Conducive to 
Approximation 

The field equations form a complete system if V{'ip) 
is known, but what we wish to do is to learn about V. 
So we follow the running of a different grouping of vari- 
ables, namely the differential equations for the set of pa- 
rameters y/e^ and ^/eyH^, with a third equation for the 
running of fl^, Eq. ([T3|) : 



--3{y/7^-^){l~e^/3), (14) 

= (e^ - em)(l - ^4>) + 27\/e0\/ev^0(-15) 



din a 



din •y/eyH 







dln( 



The reason for choosing Eq. ([T5|) for ^yevTl^ 
the simpler running equation in y/ev 

1 



-d In y/e^/d In a = 7 V^^ V^v^ 



rather than 



(16) 



is because we must allow for the possibility that the po- 
tential could be steep in the early universe, so ey may 
be large, even though the allowed steepness now is con- 
strained. For tracking models where a high redshift at- 
tractor with constant w^ exists, eyil^ tends to a con- 
stant, and is nicely bounded, allowing better approxima- 
tions. 

(The 4> ^ —(j) symmetry allows us to fix the ±.^£0 



and ±yjev7l^ ambiguity to the positive sign. For phan- 
tom energy, we use the negative sign, in Eqs. ([TH - fro)) . 
In these equations, we have restricted ourselves to fields 
that always roll down for quintessence, or up for phan- 
tom. The interesting class of oscillating quintessence 
models (jKaloper fc Sorbol [20061 : FPutta fc Scherredl2008| ) 
are thus not considered.) 

These equations are not closed, but depend upon a 
potential shape parameter 7, defined by 



I'-- 



d'^XnVjdiP- 
(9 In y/a In 1/^)2 



(17) 



It is related to the effective mass-squared in H^ units 
through 

va^^^lH'^ - 6(1 + 7)eyf^0(l - e*/3) , (18) 

Determining how 7 evolves involves yet higher derivatives 
of In y , ultimately an infinite hierarchy unless the hierar- 
chy is closed by specific forms of V . However, although 
not closed in y^evTI^, the equations are conducive for 



finding an accurate approximate solution, which we ex- 
press in terms of a new time variable. 



y — Y ^'0app ; ^ '0app 



P0oq 



(g/Qeq)^ 



Pm{a) -f /90oq 1 + (a/acq)3 ' 

(19) 
where the subscript "eq" defines variables at the redshift 
of matter-DE equality: 

^eq ^^ \PmO/ Pni,cq} 5 Pm,cq — P^cqi ^'0cq — ^^ni,cq — ^/ ■^ ■ 

(20) 
Thus y is an approximation to yfi^, pivoting about the 
expansion factor acq in the small e^ limit. By definition 

yoq = y|a=ae, = VV^- 

2.3. The Parameterized Linear and Quadratic 
Approximations 

We now show how working with these running vari- 
ables leads to a 1-parameter fitting formula, expressed in 
terms of 

e, = ±ev\a=a^, ■ (21) 

The 2-parameter fitting formula adds the asymptotic 
EOS factor 

£4100 = i£v^4>\oo ■ (22) 



The minus sign in Eq. (j21l) and (|22p is a convenient way 
to extend the parameterization to cover the phantom 
case. Here co refers to the a <C 1 limitQ The asymptotic 
equality of our 2 key variables is derived, not imposed, 
a consequence of the attractor. In addition, if the initial 
field momentum is far from the attractor, it would add 
another variable, but the damping to the attractor goes 



like 



,-6 



in Erf,, hence should be well established before we 



get to the observable regime for trajectories. 

The 3-parameter form involves, in addition to these 
two, a parameter which is a curious relative finite differ- 
ence of dyjcy^l^/dy about ycq/2: 



c 



_ dy/evTl^/dylcq - dy/ev^/dy\r 
dy/ev^/dy\cq + d^ev^^/dy\r 



(23) 



The physical content of the "running parameter" Qs is 
more complicated than that of e^ and e^oo- It is re- 
lated not only to the second derivative of In F, through 
7, thus extending the slope parameter £<, to another or- 
der, but also to the field momentum. As we mentioned in 
the introduction, the small stretch of the potential sur- 
face over which the late-inflaton moves in the observable 
range make it difficult to determine the second derivative 
of In V from the data (§ [5] In thawing models, for which 
the field momentum locally traces the slope of In V , the 
dependence of Cs on the field momentum can be elim- 
inated. However, if the field momentum is sufficiently 
small, W0 does not respond to d"^ \x\V / dcjp . We discuss 
these cases in § [6l 
What we demonstrate is that a relation linear in j/. 



yjev^<l, ~ yj£<ijoo + (V^ - V2^0oo)2/ , (24) 

'" For tracking models at high redshift, the asymptotic e<^oo 
is <x tm, and so varies by a factor of 3/4 from the radiation- 
dominated era to the matter-dominated era. In practice we ac- 
tually use |e0ool/sm as the parameter in our MCMC calculations, 

hence — 2S-emSgn(es) replaces e^ao in our w^ parametrization. 



6 



is a suitable approximation. It yields the 2-parameter 
formula for e^, whicli maintains the basic form linear in 
parameters, 



with 



/e^ ~ v^£>oo + (v^ — V 2e<)ioc ) F I 



F{x)^ 



VTT^ In [.T-^/2 ^ yr^^ 



(25) 



(26) 



2-3/2 j^3 

The 1-parameter case has ^fe^ set to zero, hence 

the simple i/eyT^ w yjejhi^ and e^ = esi^^(a/acq), 
which we regard as the logical physically-motivated im- 
provement to the conventional single-wo parameteriza- 
tion, 1-f- w = (1 + u;o)F^(a/aoq)/F2(l/aoq). 

The approximation for the 3-parameter formula adds 
a quadratic correction in y: 



yjev^(j] ~ y^£cl>co + {\/£s~ \/'2e^ac)y 



1 - Cs 1 



JL 

2/cq 



and a more complex form for the DE trajectories 

-^£0 = -\/£<^oo + (\/£l^\/2£</.oo) 



(27) 



*cq 



ana J \a. 



with 



F2{x) = V2 



In 1 



F{x) 



■cq 

(28) 



(29) 



To cover the phantom case, we put absolute values ev- 
erywhere that Es and e^oo appears, and multiply 1 -I- Wdc 
by the sign, sgn (e^) = sgn (e^oo)- 



2.4. 



Asymptotic Properties of V{(j)) & J£,^c 



The class of quintessence/phantom models where the 
field is slow-rolling in the late- universe accelerating phase 
is of primary interest. Qualitatively this implies e^ <C 1 
at low redshift, but we would like to have a parametriza- 
tion covering a larger prior space allowing for higher 
e<^ and IdV/^lnal, and let observations determine the 
allowed speed of the roll. Therefore we also include 
moderate-roll models in our parametrization, making 
sure our approximate formula covers well trajectories 
with values of |1 -I- vj^\ which extend up to 0.5 and e^ 
to 0.75 at low redshift. 

At high redshift the properties of the potential are 
poorly constrained by observations, and we have a very 
large set of possibilities to contend with. Since the useful 
data for constraining DE is at the lower redshifts, what 
we really need is just a reasonable shutofF at high red- 
shift. One way we tried in early versions of this work was 
just to cap €0 at some e^oo at a redshift well beyond the 
probe regime. Here we still utilize a cap as a parameter, 
but let physics be the guide to how it is implemented 
so that the e^ trajectories smoothly join higher to lower 
redshifts. 

The way we choose to do this here is to restrict our at- 
tention to tracking and thawing models for which the 
asymptotic forms are easily parameterized. Tracking 
models have an early universe attractor which implies 
€0 indeed becomes a constant e^oo, which is smaller than 



tra- If we apply this attractor to equation (J14p . we obtain 
•^eyil^ is constant at high redshift, equaling our ^e^oo- 
By definition, thawing models have e^oo = 0. 

Consider the two high redshift possibilities exist for the 
tracking models. One has e^ = £„, which, when com- 
bined with equation ()15p , implies the potential structure 
parameter 7 vanishes, hence one gets an an exponential 
potential as an asymptotic solu tion for V ^ a s discussed 
in Copcland ct a l. (19 98); Liddl e fc Scherreii fl999') and 
[C op eland et al.l (J2006() . Another possibility has e^ < Cm , 
hence from Eq. ([15]) we obtain 7 = (e^ — e0)/(2e0) is 
a positive constant. Solving this equation for 7 yields 
a negative power-law potential, Y ex ■0"^''''. For both 
cases, we must have an asymptotically constant 7 > to 
give a constant 



^0O 



,/(l + 27oo) 



(30) 



where 700 is the high redshift limit of the shape param- 
eter Eq. (fT7)) . This also shows why e^oo/f-m is actually 
a better parameter choice than e^oo, since the ratio is 
conserved as the matter EOS changes. 

The difference between tracking and thawing models is 
not only quantitative, but is also qualitative. For track- 
ing models, the asymptotic limit ^e^oa has a dual inter- 
pretation. One is that the shape of potential has to be 
properly chosen to have the asymptotic limit of the right 
hand side of 



^/f 0OO -^ yey^^lc 



IrflnF 

"2 dV 



^i 



(31) 



existing. We already know how to choose the potential 
~ it has to be asymptotically either an exponential or 
a negative power-law. Another interpretation directly 
relates ^£^00 to the property of the potential through 
Eq. ([50)) . In the thawing scenario, Eq. P^ is trivial, and 
the shape parameter 7 is no longer tied to the vanishing 
£000, though Eq. ([3T|) still holds. 

For phantom models the motivation for tracking solu- 
tions is questionable. We nevertheless allow for recipro- 
cal e^-trajectories as for quintessence, just hipping the 
sign to extend the phenomenology to 1 -I- w^ < 0, as has 
become conventional in DE papers. What we do not do, 
however, is try to parameterize trajectories that cross 
e^ = 0, as is done in the DETF wo-Wa phenomenology 
(see Fig. [1]). 

2.5. The Two-Parameter w^{a\es,£,f,oo)-Trajectories 

In the slow-roll limit, ey does not vary much. As 
mentioned above, we use £3 = ±e\/_eq evaluated at the 
equality of matter and dark energy to characterize the 
(average) slope of \a.V at low redshift, and ^mo or h as 
a way to encode the actual value of the potential Vb at 
zero redshift. To model y(0) at high redshift for both 
tracking models and thawing models, we use the inter- 
pretation ([5l]) characterized by the "tracking parameter" 
£000 = {^v^4>)\oo, which is bounded by the tracking con- 
dition < £000 Am < 1- 

Eq. P^ shows that to solve for e0(a) one only needs 
to know ^JeyTl^ plus an initial condition, ^e^oo, but 
that too is just y^eyH^ in the a ^ limit. We know 
as well ^Jev^^, at Oeq, namely a/es/S. If the rolling is 



quite slow, ey will be nearly e^, and fid, will be nearly 



n 



0app — y 



^, hence 



Vey^ ~ V^s^</,app = \/£7y , (32) 

the 1-parameter approximation. But we are also assum- 
ing we know the y — boundary condition, ^/ev^l^\oo 
as well as this ycq value. If we make the simplest linear-y 
relation through the two points, we get our first order 
approximation, Eq. (p4)) . for y^eyri^. 
To get the DE EOS w^, we need to integrate Eq. p4)) . 



with our -i/eyTI^ approximation. In facilitate this, we 
make another approximation. 



{^/evTXp - y^){l - £0/3) « (^£1/^0 - ^/e^) , (33) 
which is always a good one: at high redshift the track- 



ing behaviour enforces ^Jey^tp — sf^ — ^ 0; and at low 
redshift £0/3 <C 1. The analytic solution for ^/e^ retains 
the form linear in the two parameters, yielding Eq. (|25p . 
and hence the w^ approximation 



O Q-cq 



(34) 



where F is defined in 

The three DE parameters aeq, e^cx), and £s are related 
to 57„iO through the constraint equation 



exp 




— = 1 - f^. 



(35) 



ter from the trajectories. The zeroth order solution for 



a, 



only depends upon fimo- 

1/3 
cq ~ 



l-f7 



mO 



(36) 



In conjunction with the approximate 2-parameter w^, 
(j34p . this Ceq completes our first approximation for DE 
dynamics. 

2.6. T/ie Three-parameter Formula 



The linear approximation of y^ev^<piy) and the ze- 
roth order approximation for acq rely on the slow-roll as- 
sumption. For moderate-roll models (jH-w^l > 0.2), the 
two-parameter approximation is often not sufficiently ac- 
curate, with errors sometimes larger than 0.01. We now 
turn to the improved 3-parameter fit to w^; we need to 
considerably refine acq as well to obtain the desired high 
accuracy. 

The quadratic expansion Eq. ([27]) ofy^fyll^ in y leads 
to Eq. (PSI) for ^£0, in terms of the two functions F and 
F2 of a/oeq. Since Cs term is a "correction term" , we im- 
pose a measure restriction on its uniform prior by requir- 
ing ICsl < 1- Thus W0(a; Ooq, Es, e^oo, Cs) follows, with the 
phantom paths covered by £0,phantom(a;es,e0oo, Cs) = 

Sgn t^£s j£0jquintessenceV^; ps | ; |^0oo | ; ^s )• 



We also need to improve Oeq, using the constraint 
Eq. (l35t . We do not actually need the exact solution, 
but just a good approximation that works for flmo ^ 0.3. 
For example, the following fitting formula is sufficiently 
good (error < 0.01) for 0.1 < flmo < 0.5: 

, .0 /■ \ ln[r2™o/(l-^mo)] ,o7^ 

3-sgn(es)d 

where the correction to the index is 



S = 



le^ool + [0.91 - 0.78f^™o + (0.24 - 0.76f7™o)Cs 
2 



|e0oo| + (O.53-O.O9Cs) 



2|e0oo| (38) 



3. EXACT DE PATHS FOR VARIOUS POTENTIALS 
COMPARED WITH OUR APPROXIMATE PATHS 



obtained by integrating Eq. (I13p from a ~ Ocq, where 
by definition il^^cq = 1/2, to today, a = 1, hence 
flmo is actually quite a complex parameter, involving 
the entire 100(0; aoq,es,£0oo)/a history, as of course is 

^{ftmO,£s,£<poo), which we treat as a derived parame- are 



Equations (|28|) and (1371138^ define a three-parameter 
ansatz for w^ = — 1 -I- 2/3£0 (with an implicit fourth pa- 
rameter for the energy scale, flmo)- We numerically solve 
w^(a) for a wide variety of quintessence and phantom 
models, and show our 'w^{a) formula follows the exact 
trajectories very well. This means we can compress this 
large class of theories into these few parameters. 

The dark energy parameters e^, e^oo and Qs are cal- 
culated using definitions ([21]) . ([22]) and ([23]) . We choose 
z = 50 to calculate the high-redshift quantities such as 
syU-ci- Unless otherwise specified, the initial condition 
is always chosen to be at In a = —20, roughly at the time 
of BBN, at which point the initial field momentum is set 
to be zero (although it quickly relaxes). The figures ex- 
press 4) in units of the reduced Planck Mass Mp, hence 
\/2ip. The initial field value is denoted as 0ini 



In the upper panel of Figure [5] we show that once the 
memory of the initial input field momentum is lost, our 
parameterization fits the numerical solution for a nega- 
tive power-law potential quite well over a vast number of 
e-foldings, even if 1 -|- ui is not small at low-redshift. 

We know that in order to achieve both w^ ~ Wm at 
high redshift, which can alleviate the dark energy coin- 
cidence problem, and also slow-roll at low redshift, the 
negative power-law potential (or exponential potential) 
needs modification to fit the low redshift. One of the best 
known examples of a potential that doe s this is motivated 
by supergravity ([Brax fc Martin|[l999[ ): 



V{<i^) = V, [^ 



exp 



2A//n 



(39) 



where a > 11. An example for a SUGRA model is shown 
in the middle panel of Figure [5] it fits quite well for the 
redshifts over which we have data, with some deviation 
once £0 exceeds unity at high redshift. 

Another popular tracking model is the double expo- 
nential model, an example of which is given in the lower 
panel of Figure [H 

In Figures [31 we show how robust our parametrization 
is for slow-roll and moderate-roll cases, taking examples 
from a variety of examples of popular thawing models. 
The horizontal axis is now chosen to be linear in a, since 
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Fig. 2. — Examples of tracking models. The solid red lines and 
dot-dashed green lines are numerical solutions of wj, with different 



if>ini- Upper panel: 
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and AniiEr. 
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blue lines are w{a) trajectories calculated with the three-parameter 
w{a) ansatz 1281 1. The rapid rise at the beginning is the a~^ race of 
ej, from our start of zero initial momentum towards the attractor. 



there is no interesting early universe dynamics in thawing 
models. The different realizations of the w^ trajectories 
are produced by choosing different values for the poten- 
tial parameters - A for the upper panel, and Vq for the 
middle and lower panels. The initial values of are cho- 
sen to ensure ft^{a = 1) = 1 — J7,„o is satisfied. We see 
that in general our w4>{a) parametrization works well up 
to |1 -|- w| ~ 0.5. The case shown in the bottom panel is 
the potential of a pseudo Nambu Goldstone boson, which 
has been much discussed for early universe inflation due 
to (f> being an angular variable with a 27r shift symmetry 
that is easier to protect from acquiring large mass terms, 



Fig. 3. — Examples of thawing models. The solid red lines are 
numerical solutions of lo^. The dashed blue lines are calculated 
using the three-parameter to (a) formula H28| l. See the text for more 
details. 



and has been invoked for lat e universe inflation as well 
fe-g.- IKaloner fc Sorbol (l200fih l. 

In Figure m we demonstrate how the parameter (^ im- 
proves our parametrization to sub-percent level. In the 
slow-roll regime, the Cs correction is small, and both the 
two- and three-parameter formulas can fit the numerical 
solution well. But the ^^ correction becomes important 
in the moderate-roll regime. 

The last example, shown in Figure [H is a phantom 
model. 

4. OBSERVATIONAL CONSTRAINTS 

In this section we compile the updated cosmological 
data sets and use them to constrain the quintessence and 
phantom models. 
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Fig. 4. — An example showing how the f^ parameter improves 
our parametrization in the moderate-roll case. The solid red lines 
are numerical solutions of w^. The dashed blue lines are ui(a) 
trajectories calculated with the three-parameter ansatz l|28|l . The 
dotted green lines are two-parameter approximations obtained by 
forcing ^s = 0. 
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Fig. 5. — An example of phantom model. The solid red lines 
are numerical solutions of w^. The dashed blue lines are w{a) 
trajectories calculated with the three-parameter w{a) formula ||28|I . 

4.1. Current Data Sets Used 

For each of the data sets used in this paper we either 
wrote a new module to calculate the likelihood or modi- 
fied the CosmoMC likelihood code to include dynamic w 
models. 

Cosmic Microwave Background (CMB) 

Our complete CMB data s ets include WMAP 
7yr (iKomatsu et~all 120101: iJarosik et all I2Q10D 



BI 



ACT (iThe ACT Collaboration et all [2OIOI) 
CEP (IChiang et alJ I2010D. QUaD (ICastro et al 
2009|) . ACBAR (jR eichardt et al." 2009; Kuo et al 



2007t IRunvan et al. 2003; Goldstein et al. 2_00 
CBI (ISiev ers et al. 2007; Rcadhead et al. 2004 
Pearson et al.,,200a). B OOME RANG (Jones et al 
20061: iPiacentini et all 2006: M ontrov et al.ll2006[ ) 
ySA_jpickinson etalj i2004il . and MAXIMA 
(|Hananv et al., . 200"0l ). 



For high resolution CMB experiments we need 
to take account of other sources of power be- 
yond the primary CMB. We always include the 
CMB lensing contribution even though with cur- 
rent data its influence is not yet strongly detected 
(jRei chardt et al. 2009) . For most high resolution 
data sets, radio sources have been subtracted and 
residual contributions have beeii marg inalized over 
(jThe ACT Collaboration"eraII [20Ta ). However, 
other frequency-dependent sources will be lurking 
in the data, and they too should be marginal- 
ized over. We follow a well-worn path for treat- 
ing the ther mal Sunyaev-Zeldovich (SZ ) secon dary 
anisotropy (jSunvaev fc ZeldovichI 119721 Il98(il ): we 
use a power template from a gasdynamical cosmo- 
logical simula tion (with heating only due to shocks) 
described in (IBond et al.l [200a IBond et alJ [20051) 
with an overall amplitude multiplier Asz, as was 
used in the CBI and ACBAR papers. This tem- 
plate does not differ by m uch from the more so - 
phisticated ones obtained bv lBattaglia et al.l ()201Clf ) 
that include cooling and feedback. Other SZ tem- 
plate choices w ith differe nt shapes and ampl itudes 
()Komatsu fc S cliak 2(3Q1 ISehgalet al.][20m ) have 
been made by WMAP and ACT and SPT, and 
as a robustness alternative by CBI. A key result 
is that as long as the "nuisance parameter" Asz 
is marginalized, other cosmological parameters do 
not vary much as the SZ template changes. We 
make no other use of Agz here, although it can 
encode the tension between the cosmology derived 
from the CMB primary anisotropy and the predic- 
tions for the SZ signal in that cosmology. 

In addition to thermal SZ, there are a number of 
other sources: kinetic SZ with a power spectrum 
whos e shape roughly loo ks like the thermal SZ 
one (iBattaglia et al.ll2010[ ): sub- mm dusty galaxy 
sources, which have clustering contributions in ad- 
dition to Poisson fluctuations. Thermal SZ is a 
small effect at WMAP and Boomerang resolution. 
At CBFs 30 GHz, kinetic SZ is small compared 
with thermal SZ, but it is competitive at ACBAR, 
QUaD's and ACT's 150 GHz and of course domi- 
nates at the ^220 GHz thermal-SZ null, but the 
data is such that little would be added by sep- 
arately modelling the frequency dependence and 
shape difference of it, so de facto it is bundled into 
the generic Asz of the frequency-scaled thermal-SZ 
template. We chose not to include the SPT data 
in our treatment because the influence of the sub- 
mm dusty galaxy sources should be simultaneously 
modelled, and where its band-powers lie (beyond 
t ^ 2000), the primary SZ power is small. 

Type la Supernova (SN) 

We use the iKessler et al.l ()2009D data set, which 
combines the SDSS-II SN samples with the data 
from the ESSENCE proj ect (Miknaitis et al.ll2Ci07l: 
IWood-Vasev et al] 120071). the Supe rnova Legacy 
Survev (SNLS) (lAstier et all 12 0061) . the H ubble 
Space Telescope (HST) (iGarn avich c t^ 119981: 
iKnop et al.l [2003t iRiess et al.l 12004 |20d^ and a 
compilation of nearby SN measurements. Two light 
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TABLE 1 

Weak Lensing Data S 


ETS 


Data sets 


Acff 


{deg2) 


n^ff (arcmin^-^) 


COSMOS'' 

CFHTLS-wide'^ 

GaBODS'^ 

RCSd 

VIRMOS-DESCART<= 


1.6 
22 
13 
53 
8.5 




40 
12 
12.5 

8 
15 



'^ IMassev ct al. (2007); Lcssourgucs ct al. (2007). 
^ IHoeks tra ct al. (2006); Schimd ct al. (2007). 
■^ IHockstra ct al. (2002): Hockstra ct al. (2002). 
'^ iHockstra ct al.. l!2002.):.Hockstra ct al.. (.2002.'). 
" IVan Waerbeke et all (120051) : [Schimd et al.l ( 120071) . 

curve fitting methods. ML CS2K2 and SALT-II, are 
employed in lKessler et al.l ((200a) . The different as- 
sumptions about the nature of the SN color varia- 
tion lead to a significant apparent discrepancy in w. 
For dcfinitcness, we choose the SALT-II-fit results 
in this paper, but caution that if we had tried to 
assign a systematic error to account for the differ- 
ences in the two methods the error bars we obtain 
would open up, and the mean values for l-\-w may 
not centre as much around zero as we find. 

Large Scale Structure (LSS) 

The LSS data are for the power spectrum of 
the Sloan Digital Sky Survey Data Release 7 
(SDS S-DR7) Lumino us Red Galaxy (LRG) sam- 
ples (Rcid ct al. 2009). We have modified the orig- 
inal LSS likelihood module to make it compatible 
with time- varying w models. 

Weak Lensing (WL) 

Five WL data sets are used in this paper. The ef- 
fective survey area A^s and galaxy number density 
Hoff of each survey are listed in Table [T] 

For the COSMOS data, we use the Cos- 
moMC plug-in w r itten by Julien Lesgourgues 
(jLesgourgues et all I2007D with our modifications 
for dynamic dark energy models. 

For the other four weak lensing data sets we use 
the co variance matrices given by iBeniamin et al.l 
(|2007D . To calculate the likelihood we wrote a Cos- 
moMC plug-in code. We take the best fit param- 
eters a,/3, zq for n{z) ex (z/zq)" exp [— (z/zq)''] , 
and marginalize over zq, assuming a Gaussian 
prior with a width such that the mean redshift 
Zm has an uncertainty of 0.03(1 -I- Zm)- We have 
checked that further marginalizing over other n{z) 
paramete rs (a and /?) has no significant impact 
(DcLope Amigo et al]l2009| ). 

Lyman-a Forest (Lya) 

Two Lyg data sets are applied: i) the set from 
iVicl ct al.' (2004") consisting of the LUQAS sample 
(jKim ct al. 2004) and the data given iuLCroft ct al. 
(I2002D: n) the S DSS Lya data presented 



[McDonald et al.l (|2005D and iMcDonald et al.l 
()2006l ). To calculate the likelihood we interpo- 
late the x^ table in a three dimensional parameter 
space, where the three parameters are amplitude. 



index, and the running of linear CDM power spec- 
trum at pivot k = 0.9h Mpc^^. 

Other Constraints 

We have also used in CosmoMC the follow- 
ing observational constraints: i) the distance- 
ladder constraint on Hubble parameter, obtained 
by the Hubble Space T elescope (HST) Key 
Project (jRiess et al.l 120091 ): ii) constraints from 
the Big Bang Nucleosynthesis (BBN): flboh'^ = 
0.022 ± 0.002 ( Gaussian prior) and »a(z = 
10^°) < 0.2 fsee ISteigmanI [20061: lO'Meara et alJ 
2006'; 'Fe rreira fc JoyceJ 119971 iBean et al.l 120011 : 
CoDcland et al.l[2006l ) ; iii) an isotopic constraint on 
the ag e of universe lOGyr < Age < 20Gyr (see e.g. 
iDauphas 2005) . For the combined data sets, none 
of these add much. They are useful if we are looking 
at the impact of the various data sets in isolation 
on constraining our parameters. 

4.2. CosmoMC Results for Current Data 

Using our modified CosmoMC, we ran Markov Chain 
Monte Carlo (MGMG) calculations to determine the like- 
lihood of cosmological parameters, which include the 
standard six, Agz and our DE parameters, wg-Wa for 
Fig. [1] and the new ones Eg, \£,j>oo\/em and Cs for most 
of the rest. Here we use |e0oo|/em as a fundamental pa- 
rameter to eliminate the dependence of e^oa on £„ and 
the sign of Eg, whereas e^oo adjusts as one evolves from 
a relativistic to non-relativistic matter EOS. The main 
results are summarized in Table [51 The basic six are: 
^fco^^i proportional to the current physical density of 
baryons; flcoh'^j the physical cold dark matter density; 
9, the angle subtended by sound horizon at "last scatter- 
ing" of the GMB, at z; ~ 1100; InAs, with As the primor- 
dial scalar metric perturbation power evaluated at pivot 
wavenumber k = 0.002Mpc~^; Ug, the spectral index of 
primordial scalar metric perturbation; r, the reioniza- 
tion Compton depth. The measures on each of these 
variables is taken to be uniform. Derived parameters 
include z^c, the reionization redshift; Age/Gyr, the age 
of universe in units of gigayears; crgi today's amplitude 
of the linear-extrapolated matter density perturbations 
over an top-hat spherical window with radius 8/i~^Mpc. 
And, of great importance for DE, flmo', and Hq in unit of 
km s"^Mpc"^, which set the overall energy scale of late- 
infiatons. 

We begin with our versions of the familiar contour plots 
that show how the various data sets combine to limit the 
size of the allowed parameter regions. Figures [6l and [71 
shows the best-fit cosmological parameters do depend on 
the specific subset of the data that is used. We find most 
cosmological parameters are stable when we vary the 
choice of data sets. The data making the largest differ- 
ences are: the SDSS-DR7-LRG data driving flm up, and 
the Lya data driving as up. We take SN-I-CMB as a "ba- 
sis", for which nma = 0.267J:j;;^i^ and CTg = 0.817J:J5-{J2?. 

Adding LSS pushes flmo to 0.292^0^010' whether or not 
WL and Lya are added or not is not that relevant, as 
Figure [H shows. Lya pushes as to 0.844lJ^;°^g. These 
drifts of best-fit values are at the ~ 2a, level in terms 
of the cr's derived from using all of the data sets, as is 
evident visually in Figure [H 
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TABLE 2 

Cosmic Parameter Constraints: ACDM, wq-CDM, wo-Wa-CDM, es-e^oo-Cs-CDM 





ACDM 


W — Wq 


W=Wo+Wa{l - a) 


track -l-thaw 


thaw 


Current Data: CMB+LSS+WL+SNla+Lya 


riboh'^ 


U.UZZO_Q0004 


n ns'^'s+oooiH 

U.UZZO_oQQQ4 


n n99^+0-0U04 

U.UZZO_Q QQQ4 


n n99S+°°00'' 

U.UZZO_Q QQQ4 
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n -1-170+0.0024 
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n 1 1 74+00029 
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T 


0.089tHJ^ 


omotroll 
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0.96t°:°l 
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•J-yo^o.oi 


n Qfi+ooi 


0.958t°:gll 


ln(10io^«) 


3.241°:°^ 


3.24i°:°^ 


24+003 

•^•^^-0.03 


q 9/1+0.03 


q 9qQ+0.033 

o.zoa_Q 033 


Asz 


0.56t°j^ 


0.56t°:ii 


0.57t°:ii 


0.56t°:ii 


0.56trii 


^m 


0.292i^:°lJ 


0.294i°°i22 


0.293t°°}^ 


0.293troll 


9qq+ooi3 
u.zyo_o.oii 


0-8 


o.844j:°:°l^ 


n »/l 1+0-026 
U-»41_o.026 


0.847tg:°^^ 


n S/1/1+0035 
U.844_oo3g 
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Ho 


69 
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-0.98ti;lt 






Wa 
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U.UU_Q29 


£</>oo /Cm 
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Cs 
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Forecasted Data: Planck2.5yr + low-z 


-BOSS + CHIME + Euclid- WL + JDEM-SN 
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U.UZZUU_o. 00007 
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e 
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T 
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omot'oTot 


o.o9olS:°°l 
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U.9m_0.002 


0.9701°:°°^ 
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1 1 r+0.008 
0.11d_Q Q08 
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0.1id_g Qgg 


q -1 -1 r+0.009 

0.110_Q QQg 
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q 1 1 r+0.009 

0.iid_Q Qgg 


\ljYi 


o.26ot^:S]Sl 


0.261+°:°°^ 


0.260+°:°°! 


9finq+00022 
U.ZOUy_o.0022 


960'S+°°°27 
U.ZOUO_o.0024 


0-8 


0.7999tS°S}? 


0.7994lg°Si 


o.8oolS:°S! 
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0.7996l°:°Si 


Zre 


io.9tS:| 
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Ho 


72 


n+0.1 
U-0.1 
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72.0in 


71.85+°:^2^ 
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-1.00+°:°} 


-, nn+003 
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n 00^+0031 


0.008+°:°^^ 
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Tracking -I- thawing models use Wdcia\£s,£<j>oo, 
stands for "not constrained" . Hq has units km 
matter-dominated regime, e^oo = 3(1 -|- u'^oo)/2 

TABLE 3 

The marginalized 68.3%, 95.4%, and 99.7% CL 

CONSTRAINTS ON Ss UNDER DIFFERENT PRIOR 
ASSUMPTIONS. 



Prior 



Constraint 



flat prior on Eg 
flat prior on \/|e7| 
thawing prior e^oo = 
slow-roll thawing e^c 
quintessence e^ > 



0.00 



+0.18+0.39+0.72 



0.00"* 



Cs = o 



-'-0.17-0.44-0.82 

-,+0.09+0.27+0.46 

-0.07-0.28-0.52 

- _ n nn_i_o-27+o.53+o.80 

fcs — U.UUi-_Q 2g-o.61-0.98 
^ _ n ni +0-26+0.50+0.70 
£s — U.Ui_o. 28-0. 59-0. 86 
£ — 00+0-18+0-39+0-76 



Cs) of Eq. ((28| . Thawing enforces e0co = 0. n.c. 
s'^Mpc"^. 6 is in radians/100. £,„ is 3/2 in the 

Figure [7] shows how the combination of complementary 
data sets constrains Eg and QmO: the two key parameters 
that determine low-redshift observables. The label "All" 
refers to all the data sets described in this section, and 
"CMB" refers to all CMB data sets, and so forth. In 

= Figure m marginalized 2D likelihood contours are shown 
for all our DE parameters. The left panel shows that the 

- slope parameter Sg and tracking parameter e,j,oo are both 
constrained. The constraint on e^ limits how steep the 
potential could be at low redshift. The upper bound of 
£000 indicates that the field can not be rolling too fast at 
intermediate redshift z ^ \. The right panel shows that 
Cs is not constrained, and is almost uncorrelated with 
£,. This is because the current observational data favors 
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Fig. 6.— The marginalized 68.3% CL and 95.4% CL constraints 
on (Tg and QrnO for the ACDM model vary with different choices 
of data sets. For each data set a HST constraint on Hg and BBN 
constraint on Q^oh^ have been used. 
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Fig. 7. — The marginalized 68.3% CL and 95.4% CL constraints 
on £s and f2mOi using different (combinations) of data sets. This is 
the key DE plot for latc-inflaton models of energy scale, encoded 
by 1 — n^o S'lid potential gradient defining the roll-down rate, ,ye7. 



slow-roll ( 1 1 -I- ui^ I ^ 1 ) , in which case the Cs correction 
in w^ is very small. Another way to interpret this is 
that a slowly rolling field does not "feel" the curvature 
of potential. In Section |6] we will discuss the meaning 
and measurement of Cs in more detail. 

Because of the correlation between Ss and e^oo, the 
marginalized likelihood of Es depends on the prior of e^oo- 
On the other hand, the constraint on Sg will also depend 
on the prior on e^ itself. For example, we can apply 
a flat prior on d\nV/d(j)\a=a^q^, rather than a flat prior 
on the squared slope e^. Other priors we have tried are 
the "thawing prior" e^oc = 0, the "quintessence prior" 
Es > 0, and "slow-roll thawing prior" e^oo = Cs = 0. The 
results are summarized in Table [S] 

In Figure [5] we show the reconstructed trajectories of 
the dark energy EOS, Hubble parameter, distance mod- 
uli, and the growth factor D of linear perturbations. 
The Wij, information has also been compressed into a few 
bands with errors, although that is only to guide the 
eye. The off-diagonal correlation matrix elements be- 
tween bands is large, encoding the coherent nature of 



TABLE 4 
Fiducial model used in 
future data forecasts 



iLifQii ilcoh h 



0-8 



0.022 0.1128 0.72 0.8 0.97 0.09 

the trajectories. As well, the likelihood surtace in all w^- 
bands is decidedly non-Gaussian, and should be charac- 
terized for this information to be statistically useful in 
constraining models by itself. The other variables shown 
involve integrals of w^, and thus are not as sensitive to 
detailed rapid change aspects of w^, which, in any case, 
our late-inflaton models do not give. The bottom panels 
show how impressive the constraints on trajectory bun- 
dles should become with planned experiments, a subject 
to which we now turn. 

5. FUTURE DATA FORECASTS 

In this section ,we discuss the prospects for fur- 
ther constraining the parameters of the new w{a) pa- 
rameterization using a series of forthcoming or pro- 
posed cosmological observations : from the Planck satel- 
lite CMB missio n CTaubei^ '20051 : IClavel fc Tauberi 120051: 
iPlanck Science T eam 2009) ; fro i n a JD EM (Joint Dark 
Energy Mission, lAlbrecht et al.l (12006')) for Type la su- 
pernova observa tions; from a wea k lensing survey by the 
Euchd satellite (jEuclid Studv Tea m 2009); and from fu- 
ture BAO data that could be obtained by combining 
low-redshift galaxy surveys with a redshifted-21-cm sur- 
vey of moderate z ~ 2. The low-redshift galaxy sur- 
veys for BAO information can be achieved by combin- 
ing a series of ground-based galaxy observations, such 
as BOSS (.Schlcgelet al.. .20071) . For the 21-cm survey, 
we assume a 200m x 200m ground-based cylinder ra- 
dio telescope llChang et al.l 120081: iPeterson et "all 120091 : 
[Morales fc Wvithd I2009D . which is the prototype of the 
proposed experiment CHIME (Canadian Hydrogen In- 
tensity Mapping Experiment). 

5.1. The mock data sets 
5.1.1. Planck CMB simulation 

Planck 2.5 years (5 sky surveys) of multiple (CMB) 
channel data are used in the forecast, with the instrument 
characteristics for the channels used listed in Table [5] 
using Planck "Blue Book" detector sensitivities and the 
values given for the full width half maxima. 

For a nearly full-sky (we use /sky = 0.75) CMB exper- 
iment, the likeliho od C can be approxi mated with the 
following formula (jBaumann et al.ll2009() : 



21nC= Y^ (21+1)/, 



i=i„ 






sky 



In 



cr, 



/jTTrjEE _|_ /jEE(-iTT _ nnTEf-iTE 



(jTTfjEE 
TTr<EE 



+ ln 



crci 



{crf\ 



(jTTqEE _ ((JTE\2 



(40) 
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Fig. 8. — The best-fit trajectory (heavy curve) and a sample of trajectories that are within one-sigina (68.3% CL). In the upper panels 
the current data sets are used, and lower panels the forecast mock data. Prom left to right the trajectories are the dark energy equation of 
state, the Hubble parameter rescaled with H^ (1 + z)^'-^, the distance moduli with a reference ACDM model subtracted, and the growth 
factor of linear perturbation rescaled with a factor (1 + z) (normalized to be unit in the matter dominated regime). The current supernova 
data is plotted agianst the reconstructed trajectories of distance moduli. The error bars shown in the upper-left and lower-left panels are 
onc-(T uncertainties of 1 + w in bands < a < 0.25, 0.25 < a < 0.5, 0.5 < a < 0.75 and 0.75 < a < 1. 

only the dominating components in the frequency bands 
we are using, i.e., the synchrotron and dust signals, are 
considered in our simulation. The fraction of the residual 
power spectra are all assumed to be 5%. 

5.1.2. JDEM SN simulation 

For the JDEM SN simulation, we use the model 
given by the Dark En ergy Task Force (DETF) forecast 
(jAlbrecht et al.1 [20060 . with roughly 2500 spectroscopic 
supernova at 0.03 < z < 1.7 and 500 nearby samples. 
The apparent magnitude of SN is modelled as 
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Fig. 9. — Marginalized 2D likelihood contours for our three DE 
parameters derived using "ALL" current observational data. The 
inner and outer contours are 68.3% CL and 95.4% CL contours. 
(We actually show \e^ao\/^m since that is the attractor whether one 
is in the relativistic em, = 2 or non-relativistic €m = 3/2 regime.) 



TABLE 5 
Plank Instrument Characteristics 

Channel Frequency (GHz) 70 100 143 

Resolution'^ (arcmin) 14 10 7.1 

Sensitivity ^ - intensity {^iK) 8.8 4.7 4.1 
Sensitivity - polarization (i^K) 12.5 7.5 7.8 



-M 



'^°Sio(^) 



25 



-/i-^z 



,Q,2 



Full width at half maximum (FWHM) assuming Gaus- 

iii beams. 

This is for 30 months of integration. 



where l^i^ = 3 and Inmx = 2500 have been used in our 
calculation. Here, (7/ is the observed (or simulated input) 
angular power spectrum, and Ci is the theoretical power 
spectrum plus noise. 

We use the mode l descr ibed in iDunklev et al.l (|2008[ ) 
and lBaumann et all (|2009f ) to propagate the effect of po- 
larization foreground residuals into the estimated uncer- 
tainties on the cosmological parameters. For simplicity, 



^-'Z- ~ fl (5near , (41) 

where Sncai is unity for the nearby samples and zero oth- 
erwise. 

There are four nuisance parameters in this model. The 
supernova absolute magnitude is expanded as a quadratic 
function AI — jjl'^z — ^'^ z^ to account for the possible 
redshift dependence of the SN peak luminosity, where M 
is a free parameter with a flat prior over — oo < M < 
-fcx); and for ^^ and ^jP , Gaussian priors ^^'"^ = 0.00 ± 
0.03/%/2 (the "pessimistic" case in the DETF forecast) 
are applied. Finally, given that the nearby samples are 
obtained from different projects, an offset jjl^ is added to 
the nearby samples only. For fi^ we apply a Gaussian 
prior /i^ = 0.00 ±0.01. 

The intrinsic uncertainty in the supernova absolute 
magnitude is assumed to be 0.1, to which an uncer- 
tainty due to a peculiar velocity 400 km/s is quadrat- 
ically added. 

5.1.3. BAO simulation 

The "Baryon Acoustic Oscillations" (BAO) informa- 
tion can be obtained by combining a series of ground- 
based low-redshift galaxy surveys with a high-redshift 
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TABLE 6 

21-CM BAO Survey 
Specifications 



Parameter 



Specification 



O.Ol/i^Mpc"^ 

15,000 deg2 

4000 

4 years 

200 m x200ni 



shot noise 
survey area 
number of receivers 
integration time 
cylinder telescope 
antenna temperature 50K 
bias 1 

21-cin survey. We assume a fiduciai gafaxy survey with 
comoving galaxy number density 0.003/i'^Mpc~^ and sky 
coverage 20,000 deg^, which is slightly beyond, but not 
qualitatively different from the specification of SDSS- 
III BOSS (Baryon Oscillation Spectroscopic Survey) 
project. The 21-cm BAO survey using a ground-based 
cylind e r radio telescope has b een s tudied by Chang et al.l 



(|200l : lPeterson et ai.l (|2009t ) and lSeo et al.l (|2009D . The 



specifications we have used are listed in Table IHl 

For more details about the BAO forecas t techn ique, 
the reader is re ferred to ISeo fc EisensteinI (|2007f ) and 
ISeo et all (l2009h . 

5.1.4. EUCLID Weak Lensing Simulation 

We assume a weak lensing survey with the 
foUowing "Yellow Boo k" EUCLID specifications 
(|Euclid Studv Teaml[2009l) : 

/sky = 0.5 , (7int)^ = 0-35 , n = 40 galaxies/arcmin , 

(42) 
where (7int)^ i^ the intrinsic galaxy ellipticity, and n is 
the average galaxy number density being observed. At 
high £, the non-Gaussia nity of the dark m atter density 
field becomes important (jDore et al.ll2009ll . For simplic- 
ity we use a ^-cutoff £max = 2500 to avoid modelling the 
high-£ non-Gaussianity, which will provide more cosmic 
information, hence the weak lensing constraints shown in 
this paper are conservative. 
We use a fiducial galaxy distribution. 



n(z) (X ( — )^ exp 



1.5 



, with zo = 0.6. (43) 



The galaxies are divided into four tomography bins, with 
the same number of galaxies in each redshift bin. The 
uncertainty in the median redshift in each redshift bin 
is assumed to be 0.004(1 + z„i): where z„i is the median 
redshift in that bin. 

In the ideal case in which the galaxy redshift distribu- 
tion function is perfectly known, the formula for calculat- 
ing weak lensing tomography observables and co variance 
matrices can be found in e.g. iMa et akl ()2006l ). In or- 
der to propagate the uncertainties of photo-z parameters 
onto the uncertainties of the cosmological parameters, 
we ran Monte Carlo simulations to obtain the covariance 
matrices due to redshift uncertainties, which are then 
added to the ideal covariance matrices. 

5.2. Results of the Forecasts 
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Fig. 10. — Marginalized 2D likelihood contours; Using mock data 
and assuming thawing prior (e^ao = 0); The inner and outer eon- 
tours of each color correspond to 68.3% CL and 95.4% CL, respec- 
tively. See the text for more details. 



The constraints on cosmological parameters from fu- 
ture experiments are shown in Table [2] below those for 
current data. The future observations should improve 
the measurement of Eg significantly - about five times 
better than the current best constraint. There is also 
a significant improvement on the upper bound of e^oo, 
which together with the constraint on Eg can be used to 
rule out many tracking models. The running parameter 
Cs remains unconstrained. 

To compare different dark energy probes, in Figure ITOl 
we plot the marginalized constraints on r2„iO and Eg for 
different data sets. To break the degeneracy between 
dark energy parameters and other cosmological param- 
eters, we apply the Planck-only CMB constraints as a 
prior on each of the non-CMB data sets. 

The forecasted BAO, WL, and SN results, when com- 
bined with the Planck prior, are all comparable. We 
note in particular that the ground-based BAO surveys 
deliver similar measurements at a fraction of the cost of 
the space experiments, although of course much useful 
coUatoral information will come from all of the probes. 

The shrinkage in the allowed parameter space is visu- 
alized in Figure [8] through the decrease in area of the 
trajectory bundles from current to future data, showing 
trajectories of 1 -I- w^ sampled down to the one-sigma 
level. We also show the mean and standard deviation of 
1 -|- w at the center of four uniform bands in a. As men- 
tioned above, the bands are highly correlated because of 
the coherence of the trajectories, a consequence of their 
physical origin. We apply such band analyses to early (as 
well as late) universe infiation that encodes such smooth- 
ing theory priors in future papers. 

FigurelHalso plots {l+z)~^^'^{H/Ho) = ^/ptoto^/ptot.o 
trajectories, which are fiat in the high redshift matter- 
dominated regime. They are less spread out than the w 
bundle because H depends upon an integral of w. One 
of the observable quantities measured with supernovae 
is the luminosity distance. What we plot is something 
more akin to relative magnitudes of standard candles, 
51og]^o {dh / di^-.Tci) , in terms of (iL;ref, the luminosity dis- 
tance of a reference ACDM model. For luminosity dis- 
tance trajectories reconstructed with the current data 
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we choose a reference model with flmo = 0.29, which 
the SDSS LSS data drove us to; for forecasts, we chose 
^mo = 0.26, what the other current data is more com- 
patible with, for the reference. Because the supernova 
data only measure the ratio of luminosity distances, for 
each trajectory we normalize dL/di^-^cf to be unit in the 
low-redshift limit by varying Hq in the reference model. 
The error bars show n are for the current supernova data 
(jKessler et al.|[2009( ). showing compatibility with ACDM, 
and, based upon the coherence of the quintessence-based 
prior and the large error bars, little flexibility in trying 
to fit the rise and fall of the data means. 

The rightmost panels of Figure [S] shows the linear 
growth factor for dark matter fluctuations relative to the 
expansion factor, {1 + z)D{z). The normalization is such 
that it is unity in the matter-dominated regime. It should 
be determined quite precisely for the quintessence prior 
with future data. 

6. DISCUSSION AND CONCLUSIONS 

6.1. Slow-roll Thawing Models, Their One-parameter 
Approximation and the Burn-in to It 

In slow-roll thawing models, to first order w{z) only de- 
pends on two physical quantities: one determining "when 
to roll down" , quantified mostly by 1 — ^mOi ^^nd one de- 
termining "how fast to roll down" , quantified by the slope 
of the potential, i.e., e^, since e^oo = 0, and Cs ~ 0: 



1 



^^fi Islowroll, thawing 



O Ctcq 



(44) 



where F is an analytical function given by Eq. (P^ . 

l/[3-1.08(l-0™o)es] 






to 



anQ flcqislowroll, thawing ^^ \ x — Q 

sufficient accuracy. We plot F'^{x) in the left panel 
of Figure [TlJ It is zero in the small a regime and 
unity in the large a regime, and at acq is F'^{1) — 

[V2-ln(l-f 72)]^ ~ 0.284. The right panel of Fig- 
ure [H] plots the derivative dF^/dx showing \dw/da\ 
maximizes around a = acq (at redshift ~ 0.4). Three 
stages are evident in the significant time-dependence of 
dw/da: the Hubble-frozen stage at a ^ acq, where w has 
the asymptotic value —1; the thawing stage when the 
dark energy density is comparable to the matter compo- 
nent; and the future inflationary stage where dark en- 
ergy dominates (a ^ Ooq) and a future attractor with 
1 -\- w -^ 2es/3 is approached, agreeing with the well- 
known early-inflation solution 1 -\- w — 2ev/S. 

Even with such models, there is in principle another 
parameter: the initial field momentum is unlikely to be 
exactly on the e^oo = attractor. In that case, (p falls 
like a~^ until the attractor is reached, a "burn-in" phase. 
(Such burn-in phases are evident in Fig. [2] for tracking 
models.) In early work on parameterizing late-inflatons 
we added a parameter a^ to characterize when this drop 
towards the attractor occurred, and found it could be 
(marginally) constrained by the data to be small relative 
to the a-region over which the trajectories feel the DE- 
probing data. However, we decided to drop this extra 
parameter for this paper since the expectation of late- 
inflaton models is that the attractor would have been 
established long before the redshift range relevant for 
DE-probes. 
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Fig. 11. — Left panel: the function F^{x) with x = a/oeq defined 
by Eq. I I26I1 . Right Paneh the derivative dF^ /dx, showing where 
ui0 changes most quickly in thawing models, namely near acq. 
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Fig. 12. — The dependence of distance moduli ^(z) on the slope 
parameter Ss for slow-roll thawing models. The prediction of ^{z) 
from a reference ACDM model with H^o = 0.29 is substracted 
from each line. The supernova samples are binned into redshift 
bins with bin width Az = 0.3. 

The 3-parameter, 2-paramctcr and 1-parameter flts all 
give about the same result for the statistics of Ss derived 
from current data. Figure [12] illustrates why we find 
Ss ~ 0.0 ± 0.28. It shows the confrontation of banded 
supernova data on distance moduli defined by 



/.i = 25 -f 5 logio 



/ di 



VMpc 



(45) 



where d^ is the physical luminosity distance, with tra- 
jectories in /i with differing values of e^ (for the h (or 
f^mo) values as described.) The current DE constraints 
are largely determined by SN, in conjunction with the 
J7,„/i2-fixing CMB. 

6.2. A WQ-Wa Degeneracy for the Slow-roll Thawing 
Model Prior 

The parameterization (j44p can be simply related to 
the phenomenological wg-Wa parameterization by a first 
order Taylor series expansion about a redshift zero pivot: 



1 



Wo 



2 9, 1 X 

:t£.^'( — ), 

O Ucq 



and 



Wa 



2 , dF^jx) 



(46) 
(47) 



a;— l/ccq 



From formula pij) . it follows that wq and Wa should sat- 
isfy the linear relation 



1 + Wo + Wail - a*) = , 



(48) 
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where 



1- 



y 1 + "oq - Acq In 


3/2 
Qcq / 


6a?q 


/ 1 + ^1+ai^ 

in 3/2 


\ 


/ Vi+<J 



(49) 



is roughly the scale factor where the field is unfrozen. 
Over the range of 0.17 < fim < 0.5, this constraint equa- 
tion (j48l) can be well-approximated by the formula 



l + wo+Wa{ 0.264 + 



0.132 



i Lji 



. 



(50) 



Such a degeneracy line is plotted in Figure [T] Rather 
than a single wq parameter, constraining to this line 
based on a more realistic theoretical prior obviously 
makes for a more precise measurement of w. On the 
other hand, once the measurements of wq and Wa are 
both sufficiently accurate, the slow-roll thawing models 
could be falsified if they do not encompass part of the 
line. Obviously though it would be better to test devi- 
ations relative to true thawing trajectories rather than 
the ones fit by this perturbation expansion about red- 
shift zero. One way to do this is to demonstrate that 
1 + Woo is not compatible with zero, as we now discuss. 

6.3. Transforming WQ-Wa andSs-Eipoo Contours into 
Wq-Woc Contours 

The DETF WQ-Wa parameterization which is linear in 
a can be thought of as a parameterization in terms of wq 
and Woo = Wa+WQ (or edc,o and edc.oo in the DE acceler- 
ation factor language). The upper panel WQ-Wa contour 
map in Fig. 1131 shows that it is only a minor adjustment 
of the wo-Wa of Fig. [1] If we said there is a strict di- 
chotomy between trajectories that are quintessence and 
trajectories that are phantom, as has been the case in 
our treatment in this paper, only the two of the four 
quadrants of Fig. [13] would be allowed. (In the WQ-Wa 
figures, the prior looks quite curious visually, in that the 
line Wa = —(1 + Wq) through the origin has only the re- 
gion above it allowed in the 1 + wq > regime, and only 
the region below it allowed in the 1 + wq < regime.) 

To make this issue more concrete, we show in the lower 
panel of Fig. 1131 what happens for our physics- motivated 
2-parameter linear model, Eq. 1341 except it is linear for 
± yjl + if^tjel , and in the "time variable" F not in a. In 
terms of w^o and w^oo, we have ^/e^ = ^/e^ + {^/e^ - 
■y/e^o) [1 — F {a / Occ^ / F {\ / Ocq)]. The prior measure on 
l+iL'^o a-nd w^oo are uniform as in the linear-a case. Once 
the prior that only 2 quadrants are allowed is imposed 
upon the linear-a case, it looks reasonably similar to the 
lower panel. 

To illustrate that the determination of our various pa- 
rameters actually depends upon an extended redshift 
regime, one can construct window functions, that is 
redshift-filters, for the parameters. These determine how 
an error in a parameter is built up by a sum of deviations 
of the observational data from the observable constructed 
using the model for Wde, using the maximum likelihood 
formula. This shows that the window functions are quite 
extended, not concentrated around zero redshift for wq, 
not concentrated at high redshift for Woo or e^oo and not 
concentrated at Ceq for Eg- An illustration of the redshift 
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Fig. 13. — The upper panel shows the marginahzed 68.3% (in- 
ner contour) and 95.4% (outer contour) constraints on wq and 
uicxD for the conventional linear DETF parainetrization, recast as 
w = Woo + {wQ — Waa)cL, Using the current data sets described 
in § El It is a slightly tilted version of the wo-Wa version in 
Fig. [T] The demarcation lines transform to just the two axes, 
with the upper right quadrant the pure quintessence regime, and 
the lower right quadrant the pure phantom regime. If the cross- 
over of partly phantom and partly quintessence are excluded, as 
they arc for our late-inflaton treatment in this paper, the result 
looks similar to the lower panel, in which the 2-parameter £a- 
£0oo formula has been recast into 2eo/3 and 2e:oo/3, via y^o = 
v^£0oo -|- (v^ ~ \/2e0^)F(l/acq). However, as well the prior 
measure has also been transformed, from our standard uniform 
desde^cc one to a uniform deode^ao one. The quadrant exclusions 
are automatically included in our re-parameterization. 



reach of the £s parameter is Fig. [12] for the Supernova 
observable, namely a magnitude difference. The specific 
window functions are very data-dependent of course. 

6.4. Cs, the Potential Curvature and the Difficulty of 
Reconstructing V{4>) 

For slow-roll thawing quintessence models, the running 
parameter Qg can be related to the second derivative of 
In V . Using the single-parameter approximation we can 



approximately calculate d^/d In a at a 



The result 
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is: 



din a 



= \/2Afpv/^ 



' 0.533VerA/p . 



(51) 
An immediate consequence is that the amount that (j) 
rolls, at least at late time, is small compared with the 
Planck mass. In the slow-roll limit, to the zeroth order, 
Eq. ([27l) can be reformulated as 



=CsV^ 



(52) 



dN '"-"°^ ^ 2V2 

Combining Eqs. (|5HI52p . and that Mp'^d'^lnV/dct)'^ = 
—dy/2£s/d(t), we obtain 

CMp^d^lnV Cd^lnV 



where 



G«- 



^,8/ ln(lW2)' 
3 I V2 



4 dV^ 



1.005 



(53) 



(54) 



For phantom models C is negative. 

We have shown in § [5] that, for the fiducial ACDM 
model (with Ss = s^oo = 0), the running parameter C^ 
cannot be measured. For a nearly-flat potential, the field 
momentum is constrained to be very small, hence so is 
the change S(j), so the second order derivative of In V is 
not probed. To demonstrate what it takes to probe Cs, 
we ran simulations of fiducial models with varying £5 
(0, 0.25, 0.5, and 0.75). The resulting forecasts for the 
constraints on e^ and C^ , are shown in Figure 1141 We 
conclude that unless the true model has a large Es ^ 0.5, 
which is disfavored by current observations at nearly the 
2(7 level, we will not be able to measure (s- Thus if the 
second derivative of In V is of order unity or less over 
the observable range, as it has been engineered to be 
in most quintessence models, its actual value will not 
be measurable. However, if \Mp'^d'^\nV/d(/)^\ ^ 1, the 
field would be oscillating. Though in their own right, 
oscillatory quintessence models are not considered in this 
work. 

6.5. Field Momentum and the Tracking Parameter e^oo 

The late universe acceleration requires the field to be in 
a slow- roll or at most a moderate- roll at low redshift. In 
tracking models, the high redshift e^ constancy implies 
the kinetic energy density follows the potential energy 
density of the scalar field: n|/2F -^ £<t,oo/{3 — e^oo)- 
The field could be fast-rolling in the early universe, with 
large e^oo and a steep potential, and indeed is the case 
for many tracking models. If the potential is very flat, 
field momenta fall as a~^, The rate of Hubble damping in 
the flat potential limit is (/> ex a~^. However, as proposed 
in most tracking models, the potential is steep at high 
redshift, and gradually turns flat at low redshift. The ac- 
tual damping rate itself is related to the fleld momentum, 
which again relies on how steep potential is at high red- 
shift. The complicated self-regulated damping behavior 
of field momentum is encoded in the tracking parame- 
ter £000 in Eq. ([28)) (see Figure [15] where the damping of 








Fig. 14.— The 68.3% CL (inner contours) and 95.4% (outer 
contours) CL constraints on Ss and (^s, using forecasted CMB, WL, 
BAO and SN data. The thawing prior {e^ac = 0) has been used to 
break the degeneracy between Ss and e^^o ■ The four input fiducial 
models labeled with red points have ijs = and, from bottom 
to top, Ej, = 0, 0.25, 0.5, and 0.75, respectively. Only for large 
gradients can ^^ be measured and a reasonable stab at potential 
reconstruction be made 




Fig. 15. — The dependence of 'W^{z) on e^oo and C^s. We have 
fixed QmO = 0.27 and Es = 0.5. The red solid lines correspond 
to <^a =0, green dashed lines (^s = —0.5, and blue dotted lines 
fs = 0.5. For each fixed (^s, the lines from bottom to top correspond 
to e^ao = 0, 0.3, 0.6, and 0.9. 

field momentum is visualized in w{z) space). At high red- 
shift the approximation (j) ^ a^'^ completely fails, but at 
low redshift the damping rate asymptotically approaches 

Since dark energy is subdominant at high redshift, the 
observational probes there are not very constraining. In 
the 3-parameter approximation, instead of an asymptotic 
£000, we could use a moderate redshift pivot Zpivot ~ 1, 
with variables e^pivot, £s, and Cs, and the asymptotic 
regime now a controlled extrapolation. Indeed, using the 
method of § 16.31 we could define everything in terms of 
pivots of e^; for example, at z = as we used in § 16. 3[ 
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and with a third about half way in between and 1, 
such as at Oeq, which is about 0.7, so Zcq ~ 0.4. Since 
ah would be functions of our original 3 parameters and 
Oeq through the specific evaluations of F and F2 at the 
pivots, the formula does not change, but the measure 
(prior on the parameters) would be different, namely to 
one uniform in the three e^pivotS. 

Whatever the specific parameters chosen, in all cases 
the lack of high redshift constraining power means that 
our reconstructed w^{a) trajectories such as in Fig. [5] 
should be regarded as low-redshift observational con- 
straints extrapolated through theory-imposed priors on 
the smooth change in quintessence paths at high-redshift. 
The tracking assumption that we have used in this pa- 
per provides one such extrapolation, but others may be 
considered. For these, we need to define a prior on 
?!'pivot and information on the damping rate of (/>, and 
this could have an impact on the marginalized poste- 
rior likelihood of Eg- We have presented an example of 
this, showing that whereas the flat prior < je^ool < Cm 
gives Es = O.OOIq;!^^, the e^oo = "thawing prior" gives 

Es = — O.OOIq 29- We note though that because DE is so 
subdominant at higher redshift, observables such as the 
Hubble parameter H{z), the luminosity distance dL{z), 
and the linear perturbation growth rate D{z) are insen- 
sitive to the high-z extrapolation of w^{z), as we showed 
in Figure |8l Thus, even for the SUGRA model shown 
in the middle panel of Figure [2] with a time- varying at- 
tractor which is not fit as well at high z, what the data 
says about it is confined to low z, and hence the results 
obtained with our formula apply to such models as well; 
in particular the late time potential steepness must be 
controlled according to the Eg constraint for them to be 
viable. 



6.6. Using Our w^ Parametrization 

Suppose that one has in mind a specific quintessence 
potential V{(l)). How can our constraints, in particular on 
e^ , be used to test viability? The first thing is to see what 
ev looks like at a function of 0. One does not know (j)eq 
of course, but it had better lie in the range for e^ allowed 
by the data. The model would be ruled out if no (j) has 
ev penetrating the allowed region. Otherwise, the field 
equation of motion has to be solved to get 0cq, and hence 
e^ from its definition. Such a calculation is inevitable 
because quintessence models do not solve the fine-tuning 
problem. For tracking models, the low-redshift dynamics 
is usually designed (i.e., are fine-tuned) to deviate from 
an attractor. For thawing models, the initial value of 
(f>oo (to which (j) is frozen at high redshift) needs to be 
fine-tuned as well. 

Although it seems unlikely that the community will 
abandon the DETF-sanctioned WQ-Wa linear a model 
and its Figure of Merit as a way to state results for 
current and future DE experiments, our easy-to-use Es 
constraint with its clear physics-based meaning could be 
quoted as well. The 2-parameter case goes well beyond 
linear-a with the same number of parameters. Because 
we learn little about the third parameter one might 
think the 3-parameter case has one parameter too many 
to bother with, but it is not really redundant since it 
helps to accurately pave the Wc/, paths. To facilitate 
general-purpose use of this parametrization, we wrote 
a fortran90 module that calculates wsialEs .EcbocXs, ^mo) 
(http://www.cita.utoronto.ca/^zqhuang/work/wphi.f90). 
The 190 code is well-wrapped so that implementing it in 
COSMOMC is similar to what has to be done to install 
the WQ-Wa parameterization to deliver more physically 
meaningful MCMC results. 



We are grateful for useful discussions with Patrick 
MacDonald, Andrei Frolov, and Dmitry Pogosyan. 



APPENDIX 
COMPARISON WITH OTHER PARAMETERIZATIONS 



The constant w model and the linear w models, w = wq + Wa{i — a) 01 w = wq + wiz with a cutoff (jUpadhve et al.l 
|2005; Linder & Huterer 2005), are the simplest w^c parameterizations widely used in the literature. An advantage is 
that they can fit many dark energy models, including those beyond scalar field models, at low redshifts. A disadvantage 
is that they fail at z > 1 for most physical models, and in this high-precision cosmology era we cannot ignore adequate 
inclusion of the z ~ 1 information. For this reason there are various three-parameter or four-parameter approximations 
proposed to improve w at higher redshifts. Some examples are simple extensions of wo-Wa' 



An expansion quadratic in (1 — a) with a Wb parameter added to wo,Wa (iSeliak et al.ll2005l : lUpadhve et al.l[2005l : 
ILinder fc HuterCTl[2005h . 

w{a) = Wo + Wa{l — a) + Wb{l — af . (Al) 



Replacement ol (1 — a) by (1 — a^), with b a parameter added to wo,Wa (jLinder fc Huterei1l2005| ). 

w{a) — Wo + Wa{l — a ) . 



(A2) 



Four parameter models include: 



w{a) — wqWi- 



r,P' 



WiaP + Woa% 
where wo,wi,as,p are the parameters (jHannestad fc Mortsell|[2003 ): 



(A3) 
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I ^ pa,/A I _ g(l-a)/A 

where wo,Wm,ac,^ are constants (jCorasaniti fc Copelandl[200^ : 

Aw 

where Wf, Aw, at,T are the parameters (jLinder fc Hutereill2005[ ). 

Many of these four-parameter models describe a w^-transition characterized by an initial w, today's w, a transition 
redshift, and a duration of transition. For the tracking models they are meant to describe, such a phenomenol- 
ogy introduces three parameters to describe essentially one degree of freedom (the tracking attractor). Instead our 
p aramet r izatio n consistently solves the tracking 'w^{a). 

iChibal (120100 has done similar work for the negative power-law tracking models, though did not model the low-redshift 
slow-roll regime, making his parameterization much less useful than our s for compari ng with data. 

Slow -roll quintessence models have been studied by, e.g., ^Crittende n et al.l (J2007D : ISchcrrcr & SenI ()2008D : IChibal 
()2009f ). but our work makes significant improvements. Firstly, we have chosen a more optimal pivot a = acq to expand 
\b.V{4i), whereas the other works expand V{4>) at either a ^ 1 or a = 1. As a result, we approximate w^ at the sub- 
percent level in the slow-roll regime |1 + w| < 0.2 compared with the other approximations failing at |1 -I- w^l > 0.1. 
Secondly, with the field momentum from the tracking regime incorporated, our three-parameter w^{a) ansatz can also 
fit the more extreme moderate-roll and tracking models with an accuracy of a few percent. Since current data allow 
a relatively large region of parameter space (|1 -I- w^\ < 0.2 at 99.7% CL), a slow-roll prior with |1 -I- w^j ^ 1 should 
not be imposed as other papers have done. Thirdly, we have extended the parameterization to cover phantom models, 
though that is arguably not a virtue. 

Pdo(-z) and H{z) reconstruction 

One has to integrate 'w^{z) once to get the the dark energy density p^{z) and the Hubble parameter H{z). Ob- 
servational data are often directly related to either H(z) or Pdeiz), and not se nsitive to Wdciz) variations. Thus 
one may prefer to direc tly parameterize pdeiz) (jWang fc TegmarMl2004l : JAlam et al. 2004) or H{z) (jSahni et al.ll2003l : 
lAlam et al.ll200 4'. 20071 on phenomenological grounds. A semi-blind expansion of H{z) or pdc{z), e.g., in polynomials 
or in bands, differs from one in Wdciz), since the prior measures on the coefficients are radically altered. By contrast, 
a model such as ours for Wdc{z;Ss,s^ooXs) characterized by physical parameters is also a model for H{z;es,Sci,oc,Cs) 
or pdc{z;es,e^ooXs), obtained by integration. 

Where the freedom does lie is in the prior measures imposed on the parameters, a few examples of which were 
discussed in Section |H 

V{(j)) reconstruction 

For thawing models, our parameters Ss and (s are based on a local expansion of In 1/(0) at low redshift. One should 
obtain simila r resul ts by directly reconstructing the local V {(/)). An attempt at local V{(p) reconstruction was done in 
iSahlen et al.l ()2005[) . where a polynomial expansion was used, 

V{c^)^Vo + Vi^ + V2^^ + ... , (A6) 

with (/> chosen to be zero at present. A result highlighted in that paper is that there is a strong degeneracy between 
Vi and (pf). We shall see that this degeneracy is actually a result of the prior, and not something that observational 
data is telling us, as we now show. 

For simplicity we assume that flmo and h are known, hence the parameter Vq is hence fixed. Although Vi is defined 
at the pivot a = 1, it is approximately proportional to y/T^ since ey varies slowly. Because 4>o is a function of l + w\z=o, 
the degeneracy between Vi and (po is roughly the degeneracy between Ss and e^g, an expression of the low-redshift 
d ynamics being ma i nly de pendent upon the slope of In^. This is explicitly shown in Figure ITS] 

iHuterer fc PeirisI ()2007D tried to do V{(j)) reconstruction using what was then recent observational data. They 
assumed zero field momentum at high redshift (zstart = 3 is used), which is essentially the thawing prior, and used RG 
flow parameters ey, rjy, etc. expanded about the initial redshift ^gtart = 3. The constraint they find on the parameter 
ev|z=3 is very weak. From our work, this is understandable because the high- redshift dynamics are determined by the 
tracking parameter rather than by ey. They also calculated the posterior probability of ey at redshift zero, which is 
much better constrained. Their result is consistent with what we have obtained for e^. Also the large uncertainty in 
Vv\z=o the y find is similar to our result for (s- 

Recentlv iFernandez-Martinez fc Verdd (|2008D have p roposed Chebyshev expansions of V{z) or Wff,{z), a method we 



have applied to treat early universe inflation (jBond et al.i i201ft). but as we have discussed in the introduction there is 
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de facto much less inforniation in the ~ one e-folding in a that DE probes cover than the 
CMB+LSS power spectrum analyses cover. 



10 e-foldings in a that 
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